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Abstract 

This preprint contains the English translation of the paper "Relaxation dynamics of a 
quantum chain of harmonic oscillators", which was published by the author in 1980 in the 
Ukrainian Physical Journal. Comments describing its motivations, background ideas and 
results are presented as well. This paper addressed to the problem of approach to the thermal 
equilibrium in an isolated macroscopic quantum system, which was studied on the example of 
the quantum chain of weakly interacting harmonic oscillators. In the initial state, macroscopic 
energy was supplied to one oscillator (atom) in the chain. Subsequent evolution of the quantum 
state of the whole chain was determined due to the model integrability. The main subject of 
interest was the time evolution of the reduced density operators characterizing the quantum 
state of a particular atom. On the short time-scale, the energy perturbation expands along 
the chain with the velocity of the fastest phonon mode. On the long-time scale, the single- 
atom density operators display strong fluctuations around the canonical Gibbs distribution. 
These fluctuations are caused by the degeneracy of the energy level differences (presence 
of resonances) in the model of coupled harmonic oscillators. After lifting this degeneracy, 
fluctuations are suppressed providing, that the reduced density operator of each atom in the 
chain becomes very close to the Gibbs distribution at almost any time moment. 

Introduction 

The claim that an isolated (closed) macroscopic system evolves to the thermal equilibrium state 
lies in the very basis of the statistical mechanics. The question about the precise meaning of 
this statement for classical and quantum systems arises for every student studying the university 
course of statistical mechanics. With me this happens in Kharkov University in 1978. Two years 
later, I published two papers [TJH] m Russian-language journals, in which relaxation of the isolated 
macroscopic quantum system to the thermal equilibrium was studied on the example of the chain 
of harmonic oscillators. It turns out, that the set of problems considered in these papers and 
underlying ideas became the subject of considerable theoretical interest in subsequent years, to 
much extent, due to experiments on quantum dynamics in the ultracold atomic gases [3J HI [5]. 
This relates, in particular, to the quantum quench problem [HI EJ [HI El EH HU U2]> an d to the 
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equilibration scenario in quantum systems based on the Eigenstate Thermalization Hypothesis 
(ETH) [131 QU and on the idea of canonical typicality Q2J QH E0 OH EH HQ] . 

While the second paper [5] was later translated into English, the first one pQ remained without 
translation. The aim of the present retrospective publication is to present the English translation 
of paper pQ, complemented with comments relating to its motivations and results. 

In the both papers [TJ [5] , the isolated linear chain of N 3> 1 weakly coupled quantum harmonic 
oscillators ('the atoms') was considered. The model Hamiltonian reads as 

jv-i 1 

H = E 2 [ p * + " 2 {qk ~ qk+i)2 + Q2 ° q2 ^ ' n <<: n °- (1) 

Here index fc = 0, . . . ,N — 1, N enumerates the atoms in the chain, and pk denote the special 
coordinate and momentum of the k-th atom, periodical boundary conditions are implied. At 
t = 0, the chain was taken in the inhomogeneous pure state |^(0)), in which the macroscopic 
energy E = DTSlo: with n = 71/ N ~ 1 was supplied to one atom with k = 0. At t > 0, the chain 
evolved into the state \ip(t)), which was determined due to exact integrability of the model (fTJ). 

In the first paper pQ, the time evolution of the reduced density operator p(k,t) of the atom 
located at the A;-th site of the chain was calculated and studied, 

p(k,t)=Tr^p(t), (2) 

where k = 0, ...,N — 1, p(t) = \ip(t))(i/j(t)\ is the projector-type density operator of the whole 
chain, and denotes the trace over the space corresponding to the rest N — 1 atoms. 

On the short-time scale, the slow spreading of the energy from the initially excited atom along 
the chain with the velocity f2 2 /f2o was found. On the other hand, the time average of the reduced 
density operators p(k,t) of each atom in the chain over the infinite time interval was found to 
approach to the Gibbs distribution in the thermodynamic limit N — ¥ oo. It turns out, however, 
that the time-fluctuations of the single-atom reduced density operators p{k, t) remain considerable 
on the long-time scale. Such atypical behavior was associated with the strong degeneracy of the 
energy level differences (E n — E n >) in model (fTJ). It was shown in [T], that lifting this degeneracy 
by arbitrary weak nonlinear interaction provides suppression of these fluctuations and guaranties 
approach of the single-atom density operators p(k, t) to the Gibbs distribution at large t. Though 
this result was proved only for the case n < 1, it is very likely, that it holds also for n > 1. 

It was shown in [T] as well, that for a huge majority of eigenstates |n) of the total Hamiltonian 
([1]), the reduced density operators p(k, n) — Tr^- 1 |n)(n| corresponding to each individual k-th. atom 
in the chain are equal to the Gibbs distribution in the thermodynamic limit. In other words, a 
typical exact pure stationary state \n) of the chain Hamiltonian would be perceived as the thermal 
equilibrium state by the observer, which can see just one atom in the chain. This result agrees 
with the ETH, which was put forward later by Deutsch [13] and Srednicki [14]. 

In paper [2] the results of [I] were partly extended to the subsystems consisting of N' neigh- 
bouring atoms in the chain, 1 < N' <C N. Starting from the same initial pure state |^(0)), the time 
evolution of the reduced density operator p^ s \t) of the subsystem containing N' atoms was calcu- 
lated. Qualitative behavior of p^ s '(t) at large t was shown to be similar to that of the single-atom 
reduced density operator p(k,t) described in pQ. It should be noted, however, that calculation of 
the time-average in [2] was not completely rigorous, but applied a certain reasonable conjecture, 
see equation (10) in [2]. 
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The following note is in order here. The quantum chain of harmonic oscillators determined by 
Hamiltonian ([lj has the complete set of N commuting integrals of motion {ni}fL 1 which charac- 
terize the number of phonons with given quasimomenta 77 = 2ttI/N, I = 1, . . . , N. Integrability of 
model (HJ) is the direct consequence of existence of these integrals of motion. The another side of 
the coin is that conservation of N integrals of motion makes impossible the complete thermalization 
of a typical initial state in this model. However: 

(a) An observer, which can see just one atom in the chain, can not distinguish phonon modes 
with different quasimomenta 77. By this reason, conservation of the phonon numbers {ni}fL x 
in the whole quantum chain does not conflict with the chain partial thermalization, which is 
registered by the observer watching at a single atom. 

(b) Anyway, complete thermalization in model ([TJ) could be still possible for some particular 
initial states of the chain. 

The rest of this preprint is organized as follows. Section [TJ describes in detail the background 
ideas on thermalization in quantum systems, which stimulated me to do the model calculations 
described in [TJ [5] . In a very brief form, these ideas were present already in the introductory parts 
of these two papers. Section [5] contains comments on the results and some calculation details 
of paper [JJ , which English translation is given in the Appendix. Section [3] contains new results 
illustrating the effect of the nonlinear interaction on the thermalization of the atom energies in the 
chain. 



1 Thermalization in isolated macroscopic quantum system 

The fundamental question about the meaning and the mechanism of thermalization in a macro- 
scopic isolated dynamical system can be divided into two parts: 

(i) What is the precise meaning (microscopic definition) of the term 'the thermal equilibrium 
state' for the isolated macroscopic system? 

(ii) In what sense the isolated macroscopic system approaches to the equilibrium state due to its 
dynamical evolution? 

Though for classical systems, at least, partial answers to these questions could be found in the 
ergodic theory, the situation with the quantum systems was much less certain. In the statistical 
mechanics textbooks (for example, see [21)). the equilibrium state of an isolated quantum system 
with the Hamiltonian H and the energy lying in the interval between E and E + A, (A -C E) 
is usually defined through the equal a priori probability postulate, leading to the microcanonical 
density operator p A (E,H), 

p A {E,H) = - 3) 

Tr g A {E, H) 

g A (E,H)=6(H-E-A)-6(H-E), (4) 

where Tr A denotes the trace of the operator A, and 0(x) is the Heaviside step function. 

However, it is not easy to answer to the question (ii) with this conventional definition of the 
thermal equilibrium state. Really, the matrix elements of the density operator p(t) of the isolated 
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system calculated in the basis of its stationary states only change their phases with time and 
do not approach to the matrix elements of the microcanonical distribution ([3]). Furthermore, 
evolution of the isolated system from some non-equilibrium pure state IV'(O)) is described by the 
unitary operator U(t) = exp(—itH), \ip{t)) = U{t) Corresponding density operator p(t) = 

U(t) p(0) £/ _1 (£) being the orthogonal projector operator at t = 0, 

p(0) = |V(o))<V(o)|, (5) 

remains, of course, to be the orthogonal projector operator at all t ^ 0, 

P (t) = \m)m)\- (6) 

Since the microcanonical density operator ([3]) is not a projector operator, the projector-type density 
operator p(t) can not approach to it in some simple sense, say, in the strong or weak operator 
topology. 

The common way used in the textbooks to avoid this problem, is to remind that the notion 
of the isolated system is a physical idealization, and in reality, all systems interact with their 
surrounding. However, though this statement is undoubtedly correct, one could hardly believe, 
that it is the true origin of thermalization, and that the latter can not be understood within the 
concept of the isolated quantum system. 

If only isolated systems are concerned, the possible way to answer question (ii) is to restrict 
the set 'physically observable' operators Aj, to some family A and to claim, that only for these 
operators Ai € A one should expect approaching of the quantum expectation values (A(t)) to the 
their (microcanonical) equilibrium values: 

(A(t)) = Tr [At ■ p(t)} ^ Tr [A, ■ p A (E, H)}, 

with subsequent fluctuations This idea was realized by von Neumann in his Quantum Ergodic 
Theorem [22j [23] . In this theorem, operators Aj were treated by von Neumann as operators of 'the 
macroscopic observables', and were assumed to commute with each other, 
[Ai, Aj] = 0, for Ai, Aj € A. 

A deep analysis of foundations of the statistical mechanics was given by Landau and Lifshitz 
in their Course of Theoretical Physics. They define the thermal equilibrium in the following way 
(see Pages 6 and 2 in [2"4"]1: 

Definition 1.1 . If a closed macroscopic system is in a state such that in any macroscopic sub- 
system the macroscopic physical quantities are to a high deqree of accuracy equal to their mean 
values, the system is said to be in a state of statistical equilibrium ( or thermodynamic or thermal 
equilibrium) . 

The subsystems here are implied to be small compared with the whole system and weakly inter- 
acting with surrounding. Such subsystems will be called as 'allowed subsystems' in the sequel. 

In the above definition, thermalization behavior is required for the reduced family of 'the 
macroscopic observables', which relate to the macroscopic subsystems of the isolated system. 

Taking definition 11.11 as a starting point, one could arrive to a nice physical scenario, in which 
thermalization in the isolated macroscopic quantum system could result solely from its quantum 

1 Notation f(t) ~» C will be used to indicate that the value of the time-dependent quantity f(t) is close (equals 
in the thermodynamic limit) to C at almost any moment t. 
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mechanical evolution. The first step was to drop two underlined words : 'macroscopic' in the above 
definition. This extended the family A of 'thermalizable' physical observable to all physical ob- 
servables, which characterize the state of allowed subsystems, either macroscopic or not. Then, 
definition 1 1 . 1 1 becomes equivalent to the following one: 

Definition 1.2 J^j. A closed macroscopic system with Hamiltonian H at the time moment t is in 
the state of equilibrium if any sufficiently small, possibly macroscopic subsystem of it that is weakly 
coupled to the surroundings has the Gibbs distribution. 

Let us denote by a some subdivision of our isolated macroscopic quantum system into the allowed 
subsystem {s, a} with the Hamiltonian H s a , and the rest of the system {r, a}, which plays the role 
of the thermal bath. The Hilbert space C of the whole system can be decomposed into the tensor 
product C — C s . a ® C r . a of the Hilbert spaces £ s ,a and £ r . a corresponding to the subsystem, and 
its surrounding, respectively. 

Definition 1.3 An isolated macroscopic quantum system will be called to be in the (partial) ther- 
mal equilibrium with respect to the allowed subsystem {s,a}, if the reduced density operators 
p(s,a) ^-j j j.j^ s su }) S y S i ern i s t a high accuracy equal to the canonical density operators p@ 
given by the Gibbs formula 

P ( a a) =Z- 1 a exp{-f3(a)H s , a ], (7) 
Z s , a =Tt^ exp[-p(a)H B , a ], 

where Tr^'"- 1 denotes the trace over the Hilbert space corresponding to the quantum states of the 
subsystem {s, a}, and /3(a) is the inverse temperature which could depend on the subsystem {s, a}. 
Note, that the reduced density operator of the subsystem p( s ' Q )(i) is related with the density 
operator p(t) of the whole system as 

p ^(t) = Tr^p(t), (8) 

where Tr <r ' Q ' ) denotes the trace over the thermal bath space C r ^ a . 

It is clear, that the macroscopic quantum system stays in the global thermal equilibrium in the 
sense of definition 11.21 if and only if it is in the partial thermal equilibrium with respect to all its 
allowed subsystems {s, a} with the same inverse temperature /?. 

Definition 1 1 . 2 1 has several important advantages. First, it does not refer to the notion of the mi- 
crocanonical ensemble and to the equal a priori probability postulate. Second, with definition 1 1.21 
the answer to the question (ii) becomes straightforward and leads to the natural treatment of the 
thcrmalization phenomenon in quantum systems. 

Proposition 1.4 Thermalization in the macroscopic isolated quantum system means, that after 
the relaxation time, the matrix elements of the reduced density operators p( s,a '(t) of all its allowed 
subsystems {s,a} in the basis of their stationary states l^rf'"') should approach (up to further 
fluctuations) to those of the canonical density operators (?[), 

( ^,a)| p (a,a) (t) | $ (.,a) ) _ ^-1 ^(.^..a)) ^ (9) 

where Z s , a = £ expH^"'), 
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and £n ' , n = 1, 2, . . . is i/ie energy spectrum of the subsystem {s, a}. 

Third, one can easily see, that the problem mentioned in the beginning of this Section does not 
arise with definition QTU Really, if the whole system is in a pure state described by some projector 
density operator (j6]), the subsystem reduced density operator (|8|) can correspond to a mixed state 
due to the quantum entanglement, and this mixed state can approach to the equilibrium Gibbs 
state 0. 

If a macroscopic isolated quantum system possesses the property of thermalization, one should 
expect, that a typical initial state 1^(0)) having macroscopically well determined total energy would 
approach in its unitary evolution to the equilibrium state, understood according to definition 11.21 
and then remain in this state almost at any time. 

It is tempting to take the initial state |^(0)) of the whole system as a single stationary state 
\n) of its Hamiltonian, H\n) — E n \n) . Then, the previous statement naturally leads to the ETH 
[THE], which can be formulated as follows: 

Proposition 1.5 In a thermalizable isolated quantum system, the huge majority of individual 
stationary states \n) of the whole system represent the thermal equilibrium state. This means, that 

Tr^|n)H =Z-iexpH8 (10) 

for all allowed subsystems {s,a}, and almost all stationary states \n). 

Though the described above heuristic scenario of thermalization looked rather reasonable and 
attractive for me in 1979, I had no idea how to put it on a more firm basis for a generic isolated 
quantum system. So, it was natural to try to test this scenario on some simple model of interacting 
particles. Such a test was partly realized in [TJ[2] for the model of the quantum chain of harmonic 
oscillators. 

It should be noted, that definition 11.21 provided the basis of the described above treatment of 
thermalization. It implies that the family A of thermalizable observables can be associated with all 
observables which characterize the states of allowed subsystems. Approach to the thermalization 
problem emphasizing the role of the canonical ensemble for the subsystem state as the main 
indication of thermal equilibrium in quantum systems0 was developed later by several groups. 
Important results along this direction were obtained by Tasaki [2D], Popescu, Short, and Winter 
[15] . Linden et al. [16], and by Goldstein et al. [T7l f!8] . 

2 Comments on the results of paper [lj 

In the normal mode representation, the quantum chain of harmonic oscillators reduces to the free 
phonon gas. As it was mentioned in the Introduction, conservation of phonon numbers {rir}^ 
makes impossible complete thermalization of a typical initial state |^(0)) in this model. However, 
the quantum measurements performed on one atom in the chain does not allow one to distinguish 
different phonon modes. So, it is natural to expect, that for model fl}, the weak versions of the 
thermalization property (proposition [O]) and of the ETH (proposition [E5J are still valid, in which 
Equations ([7]) and ([5J are required only for the subsystems {s, k} consisting of one k-th atom in the 
chain[f| Paper pQ, which English translation is given in Appendix, was dedicated to verification 
of this statement. 

2 The origin of this approach comes back to Erwin Schrodinger 1251 . see discussion in |17| . 

3 Note, that due to the requirement Q -C Qoi each atom in the chain weakly interacts with its surrounding and, 
therefore, represents an allowed subsystem by itself. 
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2.1 Eigenstate thermalization for free bosons 

The Eigenstate Thermalization Hypothesis (see proposition 1 1 . 5|) was the important part of the 
described above heuristic scenario of thermalization, which provided motivation for [TJ [5] . Note, 
that the ETH was not formulated explicitly in [TJ [2J, this has been done later by Deutsch [TJ] 
and Srednicki [14] . However, a weak form of the ETH for model ((TJ) was proved in [JJ. Namely, 
equation (|A.8[) validates the ETH (proposition II .5[) in model (jTJ) for the single-atom subsystems. 
It was conjectured also in |TJ, that, "perhaps, this result holds to some extent in other dynamical 
systems" . 

Parameter n = yi/N in (|A.8|) is related with the inverse temperature (3 by the Plank's formula 

" = cxp(/3fio)-l' (U) 

and Ul — Y^,h=i n i- I n the case h < 1, equation (|A.8p can be easily obtained from (|A.7[) which, in 
turn, follows from the results of Appendix 1 of paper [2j. 

To prove (|A.8|) . let us rewrite (|A.7[) for n = n' in the explicit form 

91 ^ r>\2( -\\P-J N I 

J=0 F=J V ; ^m F P tl l 1 1 

where m denotes the set of integer numbers {m;}^, which lie in the interval < m; < nj. Denote 



by Ni the number of nonzero integers in the set {n/}^. For fixed P > and macroscopic iVi, 
the leading contribution to the sum 

S 2 (n 1 ,...,n N ;P) = V IT- ^- (13) 

^ " (nj - miy.mv. 

2^, mi=P 1=1 

comes from configurations {mi}fL 1} in which m; takes values either 0, or 1. Contribution of all the 
rest terms is smaller by the factor N^ 1 . The same is valid also for the sum 

N , AT 

5 1 (n 1 ,...,n JV ;P)= ]T [J? ^TT^ = pi ( fl *) P U=o I^ 1 + z )"' = ( 14 ) 

(71/ — 777,7 ) ! 777/ ! P! J - J - 
mi=P i=l v ; i=l 

\(0.) p U=o(i + *)*- m 



Therefore, 



PV ' 1 v ' P!(9t-P)! 



5 2 (m, • . . ,njv; P) = 5i(n x , . . . , n N - P)[l + 0(A^r )] = (15) 
fyri *Y7 P 

' [i + (ivr 1 )] = 4 r [i + o(Ar 1 - 1 )] ! 



P!(9T-P)! L v 1 yj P! 

at fixed P and 01 ~ N\ ~ AT — > oo. Formal substitution of ([15]) into jT2j yields 

00 00 Plf-1 S i p - J s p 

Tr«|n)<n| = £ |J)(J| £ }u P _j V i 1 + °( N i 1 )} = ( 16 ) 

.7=0 P=J >' 



J=0 v ; 
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in agreement with (|A.8[) . If n < 1, the above formal manipulations can be justified, since the 
infinite scries in P in (II 6[) converges uniformly in N — > oo. 

At n > 1, the series in P in (|16p diverges. Nevertheless, it is also possible to prove (|A.8[) at 
n > 1 directly from (IA.7I) . though the proof is more complicated, than at n < 1. On the other hand, 
a much easer indirect way to confirm (| A.8|) exists, which is suitable for all positive n. Namely, one 
can check that the quantum expectation values of the operator b\ Jl b J § calculated with the both 
sides of this equation lead to the same results for any integer J\, Ji > 0. 

Proposition 2.1 Let {ai}^L 1 and {b k } k=1 are the two sets of N bosonic annihilation operators, 
which are related by the linear transformation 

N 

b k = ^A kl {N)a l , (17) 
i=i 

with the unitary N -dependent matrix A k i{N), and 

\n 1 ,...,n N ) = (n 1 l... r^!)- V Vi)" 1 ■ • ■ |0), (18) 

denotes the basis associated with the first set of the bosonic operators. For each N , let us choose a 
basis vector |ni(iV), . . . , nzv(iV)) in such a way, that 

\A k i(N)\ 2 ni(N) <C/N for all k, I G [1, N], with some fixed C > 0, (19) 

and the limits 



n{k) = lim V|A fci (A0| 2 "KA0, (20) 

TV— >oo — ' 

1=1 

exist for given k. Then 

lim (n N {N),...,n 1 {N)\bl Ml b*f 3 \n 1 (N),...,n N {N)) = M x ! [fi(fc)] Ml 5 Ml M 2 . (21) 

N— too 

Proof First, it is clear, that operator b\ Ml bff 2 with Mi ^ Mi has zero diagonal matrix elements 
in the basis (|18p . since this operator changes the number of bosons. So, it is sufficient to put 
Mi = M2 = M in (|21l) . Then, substitution of (ITTl) into the matrix element in the left-hand side of 
this equation followed by straightforward calculations yields 



(n N , m\b\ bf |m, ...,n N ) 
Ml 2 
h\ 2 ...l N l 2 



51 / \2 M ' 1 |2 l^fcZ! ■ ■ • A k i N \ 2 (n N , . . .,m\a% N . ..a[ h a^ . ..a l ^\m, . . .,n N ) 



h + ...+l N =M 



^ M , 12 n x \...n N \ 

hP...l N P lAkh --- AklNl (n 1 -l 1 )l...(n N -l N )\ 



^ 7772 j^lAkh ■ ■ ■ Aki 

h + ... + l N = M 

QM 

Ml 



dt M 



l\L ni (-t\A kl \ 2 ), (22) 



*=0 1- 



1=1 



s 



where L n (z) is the Laguerre polynomial, 
n\(-z) 1 



L n (z) = J2 



1=0 



l\ 2 (n-l)V 



The explicit form of (|22|) read as 

(n N , . . .,ni|6|.6fc|ni, . . . ,n N ) = ^ \A k i\ 2 n h 



N 



1=1 



(n N , . . . ,ni|6^ 2 6||ni, . . .,rijv) = 2 



(riAr, . . .,m\b£bl\m, . . . ,niv) = 3! 



N 



N 



U=l 
' N 

\ A kl? n l 



N 



.1=1 



-^|A fci | 4 (n i + nf), 
i=i 

3 N 

+ J2\ A ki\ 6 (2m+6n 2 +4nf) 



/=1 



m— 1 



(rajy, . . . ,m|6fc ^|m, . . . ,n N ) = Ml 



N 



1=1 



1=1 



M 



It is easy to understand, that only the first term in the right-hand side of the above relations 
survives in the thermodynamic limit N — > oo due to (|19j) and (|20j) . 



Let us associate with each pair of operators b\,b k the 'single-atom subsystems' {s, k} with the 

(s) t 

Hamiltonians nl = e k b k b k . Suppose now, that each such atom is in the thermal equilibrium state 
characterized by the fc-dependent inverse temperature 

1 + n(k) 



(3 k = e, 1 log 



n(k) 



This means, that the reduced density operator of the fc-th atom has the Gibbs form, 



pa(k) = 



exp [-p(k)h k > 

7 (s) 



Wi 



E 



[my 



-r. b 



^ [1 + m y+i j! 



b(. 



(23) 



(24) 



Then, for arbitrary natural Mi , Mi , one can easily obtain the expectation value of the product 
4 Ml& f 2 in this state 



Tr 



s,k) 



bl Ml bt h PG(k)} = M x \ [n(fc)] Ml 6 MlM2 , 



(25) 



which coincides with the right-hand side of (|2"Tj) . 

Thus, at large N, the huge majority of basis states (fT51) correspond to the partial equilibrium 
states with respect to the single-atom subsystems {s, k} in the sense of definition 1 1.31 The inverse 
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temperature f3k of the fc-th atom, generally speaking, depends on fc, and is given by equations (|23[) . 
(HQ]). 

Let us return to the chain of harmonic oscillator ([1} . For the bosonic annihilation operators b k 
and ai determined according to (IA.3I) . the matrix Aki(N) introduced in (fTTf reduces at O = to 
the simple form Am = N~ x / 2 exp (2iii l/N), and e k — Qq. Parameter n(fc) defined by (|2"0"|) reduces 
to the average number of phonons in the chain n — ^V~ 1 X)i=i n ' an d becomes fc- independent, 
together with the inverse temperature (f2"3|) . This proves (|A.8|) for all n > 0. 



2.2 Time evolution of the single-atom density operators 

Relaxation of the macroscopic chain of harmonic oscillators from a nonequilibrium nonuniform 
initial state |V>(0)) = (tft!) -1 / 2 fc^O) was studied in Section 2 of pQ. First, the time evolution of 
the chain quantum state \tp{t)) was determined in the leading order in the interaction constant 
fl — > 0. Then, representation (|A.12[) was obtained for the density operators p(k, t) of the fc-th atom 
at time t, which is asymptotically exact at Q — > 0. The full form of this expression reads as 



P (k,t) = -£\J)(J\ 



with 



j=o 



N 



j\{m-j)\ 



r(k,t) 



N 



2 J 



r(k,t) 



N 



(26) 



r(k,t) 



^exp [i{jik - LJ n t)] 



i=i 



(27) 



Remind, that 77 = 2t:1/N and w T are the phonon quasimomcntum and the energy, respectively, 

n 2 



[% sm 2 (r/2) 



1/2 



n 



■(i 



o{n 4 



(28) 



Parameter r 2 (k,t) is proportional to the energy E(k,t), which is located at the fc-th atom at 
the time moment t: 



E(k,t) = n (m\b{b k \m) = n Ti- (s ^[ P (k,t)blb k ] - yin r 2 (k,t)/N 2 



(29) 



At small enough t, the sum in (|27|) can be replaced by the integral. This allows one to express the 
energy E(k, t) normalized to the total energy in the chain in terms of the cylindric Bessel function 
Jk{z), 



E(k,t) 



2 (M) 

N 2 



dr 
27' 



t i(rk— uj T t) 



dr 
2^ 



exp 



-k 



n 



COS T 



[j k (tn 2 /n )] ■ (30) 

The plot of this normalized energy versus k and t is shown in Figure [TJ It is clear from (IA. 141) 
and Figure [l] that the energy perturbation spreads along the chain with the velocity v — 2 /rio, 
which is just the maximum group velocity of phonons 

v ph (r) = 

v = maxvphir) = v ph (ir/2). 



(n 2 /ft ) sinT + 0(0 4 ), 
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This feature of the relaxation dynamics known as 'the horizon effect' is typical in the quantum 
quench phenomenon [7]. 

2.3 Dispersion of the density operators 

Now let us turn to the dynamics on the long time scale. When time t becomes larger then the 
relaxation time t r = N/v = Nflo/H, 2 , the energy perturbation reaches all atoms in the chain, 
all exponential terms in the sum in right-hand side of (|2T[) completely dephase from each other, 
and this sum cannot be approximated by the integral, as it was done in ([30)1 . This leads to the 
qualitative change of of the chain evolution character from regular at t < t r , to the 'stochastic' 
regime at t > t r , which is illustrated in Figure^ 

For the atom reduced density operators averaged over the infinite time interval, the canonical 
Gibbs distribution (|A.17[) with the site- number independent temperature was obtained in [T]. 
However, derivation of this result was incorrect, since it was based there on the erroneous equation 
(A. 14a) in Page[25j see the footnote in this page. The fallacy in this equation came from the implicit 
and mistaken assumption, that the energy levels E n of the whole system are non-degenerate, 

E n ^E n ,, iin^ri. (31) 
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1.0 
0.8 
0.6 



0.4- 




Figure 2: The red line shows the time dependence of the quantity r 2 (0, t) /N 2 defined by (|27|) . which 
is proportional to the energy £"(0, t) localized on the 0-th atom, E(0,t)/ (yiflo) — r 2 (0,t)/N 2 . The 
blue line plots [Jo(tn 2 /£lo)] 2 . The model parameter values are taken as N = 30, Sl 2 /f2o = 0.01. 
Change of the evolution regime takes place at t ~ t r , t r H 2 /H,q = 30. 



Really, the energy levels (measured from the ground state energy) in model (JlJ are given by 

E n = (nu>) = n T u> T , (32) 

where summation in the quasimomentum r runs over the Brillouin zone, 

S = t£ • (Z mod N). Due to the the mirror symmetry cu T = oj_ t , of the phonon dispersion law 
(|28[) , this energy spectrum can be rewritten in the form 

E n = n Q uj Q + ^^[l + (-l) N ]+ £ («r+n-r)wr, (33) 

tEB 
0<r<7T 

which indicates the strong degeneracy of the energy levels E n in model (TT]), in contrast to (|3"Tj) . 
Accordingly, one gets for the time average exp[it(E n > — E n )\: 

!<^no,<^,<n t£B ^Tv+n-T.n'+nL > for even AT, 
a rr a <T< " " t AA AT ^ 

d n ,n'[[ tEB <V+n_ T ,n' +n> , for odd N, 

0<T<7T 

instead of equation (A. 14a). It turns out, however, that this correction does not change the final 
result (|A.17[) for the time averages of the reduced density operators almost for all atoms in the 
chain. 
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Proposition 2.2 In the thermodynamic limit, the time averages of the atom reduced density op- 
erators p(k,t) determined by equation ([A. 121) [or, eqivalently, by equations K26\) , (|27|) / are given by 
the canonical Gibbs distributions with two different temperatures, 



hm p(k,t) 

N— >oo 
<n/N=n 



IXo { i+my^ I J > ( J \ > for fc = N / 2 and even AT, 



Ej=o [TRp+i otherwise, 
where the time average is defined as 



1 f 1 

p(k, t) = hm - / 
T-s-oo T J 



dt p(k, t). 



(35) 



(36) 



Equation (|3"5j) corrects equation (|A.17I) for fc = 0, and k — N/2. 

Proof will be given only for the case of even N, since extension to the case of odd N is straight- 
forward. 

Let us introduce the (-^ + l)-dimensional phase space with the points <j> — {4> T }, where < 
(\> T < 2ir, and r = t(1) — 2ttI/N, with I = 0, 1, . . . , N/2. Let us define the operator-valued function 
p(k, </>) on this space which is obtained from (|A.11|) by the replacement 



t{nuj) -> (n<p) = n 4>o + + ^ (n T + n_ T )0 T , 

0<T<7T 

t{n'uj) -> (n'4>) = n' (f>o + ri^fa + « + n'_ T )</>r, 



0<T<7T 

and averaging over the phases {cj> T } of the 27r-periodical functions f(<j>), 

,2^ N/2 



(/(<£»< 



2=0 



2tt 



Then, it follows from f3"l]). that 



exp[ii(w(ri' — n))] = (exp[z(0(n' — n))])^. (37) 

Applying this equality to (|A.11[) . one can conclude, that the time average of the operator p(k, t) 
is equal to the phase average of the operator p(k, 4>). Therefore, the time average of the operator 
(l26l) can be written as 



J=0 



p(k,t) = J2\J){J\ Cm ([r(k^)/N] 2J 1 - [r{k,4>)/NY 



m-j 



(38) 



where 



r(k,(j>) 



N/2 
1=0 



wi(k) 



1, for I = 0, 

2cos[fcr(Z)], for I = 1, . 
(-l) fc , for/ = f. 



N 
2 



1, 
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Accordingly, the matrix elements of the operator (1551) in the basis | J) = (J!) 1 / 2 b^ I \0) take the 
form 



{J'\p(k,t)\J) = 6jj. PJ (k,t), 
PAr(fc, r) = (<5[r - r(fc, 0)])^ 



dr(r/Ny J l-(r/N)' 



PN(k,r), 



(39) 
(40) 
(41) 



After rescaling of the integration variable r = v N x in the right-hand side of (|40p. the large- JV 
asymptotics of this equation can be written as 



(42) 



pj(k,t) = J dx— expi-nx^^NpNik^X^^ + OiN- 1 ) 



Function pjv(fc, r) defined by (|4"Tj) is the probability density to find the total length r of the sum 
of N/2 + 1 vectors in the plane which have arbitrary directions and different lengths I = 

0, . . . , N/2. Straightforward calculations lead to the exact representation of this probability density 
in terms of the cylindrical Bessel function Jq(z) [cf. equation (|A.16[) in Appendix], 



p N (k,r) = r 



N/2 

dvvJ (vr) J [wi(kY 

1=0 



(43) 



At large N, the main contribution to the integral in the right-hand side of (|43|) comes from small 
v, since the integrand vanishes at v ^> TV -1 / 2 . So, in the thermodynamic limit we get from (|43|) : 



Np N (k,r)\ r=xVW = Nx j dvvJofvxVN^J 
exp[-x 2 /Y(k, N)] + O^ 1 ), 



exp 



Nv x 



Y(k,N) 



[l + O^- 1 )] = 



2x 



Y(k, N) 



where 



1 N/2 i2 

Y(k,N) = -J2m(k) 2 = \ 

1=0 I 



, iffc = Q,jV/2, 



(44) 



(45) 



Substitution of ([44]) and ((45]) into ((42]) yields finally to the result 
[nY(k,N)] J 



[l + nY(k,N) 



J+1 ( } {jT^ + OiN-i), iffc=l,...,#-l, W 



Pj(k,t) 



which, together with (f39l) proves (|35|) . 

Let us return now to the time evolution of the atom reduced density operators (|A.12|) at large 
t > t r . At a fixed time moment t > t r one can treat r(fc, t) as the absolute value of the sum of N 
unit vectors on the plane which have 'stochastic' phases. Accordingly, it is again natural to rescale 
the 'length of the random walk path' 



r(k,t) = VN x{k,t), 



(47) 
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with x(k, t) ~ 1 for almost any t > t r . Substitution of (|T7)) into (| A. 12|) leads in the thermodynamic 
limit N — > oo, 91 — > oo, OT/iV = ri to the Poisson distribution 

p(M) = £ |J)(J| [ " x2( ^ )]J exp [-m 2 (M)] (48) 
j=o J - 

for the atom density operators at almost any time moments t 3> t r . On the other hand, the time 
average (|36[) of the density operator p(k, t) over the infinite time interval was proved to approach 
to the canonical Gibbs distributions in the thermodynamic limit, see (|35[). It follows from (1351) 
and (|48p . that the reduced density operator p(k,t) of the fc-th atom never approach to its time 
average ((35)) at any time moment. Similarly, the quantum averages (Ak)(t) = Tr^ s ' k ^ [A^ p(k, t)] of 
the local observables relating to the fc-th atom should also strongly fluctuate in timeQ This 
means, in turn, that the time dispersion 



D,j.j>{p) = pj{k,t) PJ ,(k,t) - PJ (k,t) pj,(k,t) (49) 

of the density operator matrix elements p,j{k,t) — (J\p(k, i)\ J) remains considerable and does not 
vanish in the thermodynamic limit. Such type of evolution indicating luck of partial thermalization 
in model (TTJ was associated in [T] with the strong degeneracy of the energy level differences, see 
equation (|A.18|) . Really, the dispersion (pt9")l can be written as 

D JJ'(P) = E E ( J l Tr(r ' fe) [K>«ll] l J > (J'\Tr™ [K){n 2 \] \J') ■ 

ni ,71^ 712 ,7T<2 

(ni|V(0))(^(0)|ni)(^^(0))(^(0)|n 2 ) • 

{exp[zt(^„ 2 - E n , 2 - E ni + £„,)] - exp[zt(£„ 2 - E n > 2 )] exp[it{E n[ - £ ni )]} , (50) 
where E n are the energy levels (|32p in model (H}. The last line in ([501) has the structure 

{exp[it(E n2 - E n > 2 - E ni + # n > )] - exp[it{E ri2 - £?„/)] exp[it(E n(i - E ni )]\ = 

where 5 nin . 2 is the Kronecker delta, and 



^n\n2n' 1 n 2 



1, if n-y = 712 = Jlj = ^2, 

otherwise. 



The resonance term C nin2n > in > 2 in the right-hand side of (|51j) is nonzero, in particular, at n\ — n\ = 
n' 2 — ri2 = m, and as well at n-i — n\ = n' 2 — n\ = m due to (|A. 181) . It is precisely this term which 
is responsible for the strong time fluctuations of the density operator p(k,t) in model (TTJ). 

On the other hand, the resonance term C ni „ 2 „<„< vanishes in (|51|) . if a weak anharmonic 
interaction is applied providing for the shifted energy levels that 

„ „ , „ „ , either ni = n\ and no = nf , 

E ni - E n , ± E„ 2 - E n , unless <^ 1 1 j f 2 (52) 

I or rii = i%2 and n\ = n 2 . 



For the 0-th atom energy operator h k , the time fluctuations of its quantum expectation values E(0, t) are 
clearly seen in Figure [2] at tQ 2 /Qo > 30. 
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In this case, one should expect, that the time dispersion Djji{p) would also vanish in the ther- 
modynamic limit N — > oo,9T/-/V = n. For n < 1, this was indeed proved in pQ. The crucial role 
of the 'non-resonance condition' (|5^|) for ability of an isolated quantum system to thermalize was 
first established by von Neumann [53] , and confirmed later by many authors ^ZU\ EH HH] ■ 



3 Time dispersion of the quantum expectation value of the 
atom energy 

It is interesting to illustrate the dramatic effect of the nonlinear interaction on the time evolution 
of the quantum expectation value of the energy E(k,t) of a particular atom, 



E(k,t) = m)\hl s) \m), 



(53) 



where 

is the Hamiltonian of the k-th atom, and \tp(t)) is the quantum state of the chain at time t which 
evolves from the initial state (lA.lpp . 

In the purely harmonic chain (fTJ), the energy E{k,t) is determined by 



E(k,t) = mn { 



r 2 {k,t) 
N 2 ' 



Taking into account ff2T|) . one obtains from (f54| for the time averages: 

if k = 0, 

if k = N/2, and N is even, 
otherwise, 



lim E(k,t) 

N-HX> 




(54) 



(55) 



lim D(h[ s) ) 

N— >oo 



" ) 
2 ^2 



4^o n 



if k = N/2, and N is even, 
otherwise, 



(56) 



where 



D{h ( ° ] ) = E{k,t) 2 - \E{k,t) 



(57) 



Thus, the time dispersion D(h/! ) of the quantum expectation value of the atom Hamiltonian 
remains finite in the thermodynamic limit in the purely harmonic chain. 

On the other hand, in the chain perturbed by the extremely weak anharmonic interaction, the 
fc-th atom energy (|53|) can be determined as 



E(k,t) = Tr^[h^p(k,t)}, 



(58) 
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and the atom reduced density operator 

ml rm/Tif P! 2 (-i) JJ - J x , I 



E exp{^[fc(nr) - ^ exp{i[-fc(n'r) +tE n ,]} 

(n-l)\ ^ (n'-l)\ ' ( } 

can be obtained from (|A.11[) by replacement 

(con) -> E n = (con) + SE n , 
(am') -> = (wn') + 5E n i, 

corresponding to the small shift of the energy levels. 

If the perturbed energy levels E n are not degenerate, averaging in time of with the atom 
density operator p(k, t) given by (|59p provides 



E(k, t) = fl n, 

instead of (1551) . In turn, if the energy levels E n satisfy the non-resonance condition (f5"2"]) , the time 
dispersion (|57[) can be written as 



oo oo 



D(4 s) ) = n^£ jj'DjAp), (60) 

J=0 J'=0 

where Djj' (p) is given by (|A.19j) . Due to the equality 

(-l) p - J P! fl, if P = l, 
i ~ 

several summations can be explicitly performed in (|6ip . The result reads as 



tj [ ~ 1} r =r> 11 ^ = i ' (en 

^ J\(P-J)\ )0, otherwise, V ' 



V k ' °N^+ 2 ^ l\(n - l)\{n' - l)\ k\(n - k)\(n' - k)\ 



Er 'r=Er *x=9t-l 



where integration is performed in the complex x-plane along the circle centered at origin and going 
in the positive direction, and Io{t) is the modified Bessel function of the first kind, 

j=o J - 

In the thermodynamic limit N — ¥ oo, n — 9t/iV — Const, the integral in the right-hand side of 
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1 2 3 

Figure 3: Function W(n) determined by ffGGj) . 



n 



(|62| is determined by contributions of two saddle points at x = ±iro(n), where a;o(n) is the positive 
solution of the equation 



h{2x) n 
Iq(2x) x 



= 0. 



=x (n) 



The saddle-point asymptotics of (p^)) in the thermodynamic limit reads as 



2nN 



1/2 



[W(n)] N [1 + OiN- 1 )], 



where 



f(x)=log[I (2x)}-2nlogx, 



and 



n\ 2fi I [2cc (n)] 



(63) 



(64) 



(65) 



(66) 



The plot of function W(n) is shown in Figure [3] This function has the following asymptotic 
behavior 



W(n) 



(n/e 2 )™ [1 + n + 0(n 2 )], n^+0 : 
(47rn)- 1 / 2 [1 + C^rT 1 )], n -> +oo, 



(67) 
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and monotonically decreases between values W(0) = 1 and lim^oo W(n) = in the half-axis 
< n < oo. Therefore, it follows from (|64|) . that the time dispersion (|57p of the quantum 
expectation value of the atom energy vanishes in the thermodynamic limit N — ¥ co for all fixed 
fi > 0. 
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Appendix 



Relaxation Dynamics of a Quantum Chain 
of Harmonic Oscillators Q 

S. B. Rutkevich 
Kharkov University, Kharkov, Ukraine 

Abstract 

The aim of the present paper is to trace the appearance of some statistical prop- 
erties in an exactly solvable dynamic system - a quantum chain with weak harmonic 
interaction. The time dependence of the statistical single particle matrices is derived 
when the system relaxes from the non-equilibrium pure state. It is shown that the 
thermodynamic limit of the matrices in question is the Gibbs distribution. Moreover, 
the presence of any weak anharmonic interaction in the system causes each oscillator 
to have a canonical distribution almost at any moment. 

Introduction 

The statistical mechanics, besides the dynamical laws, is based on certain additional probability 
hypothesis (the micro-canonical distribution postulate, the correlation weakness principle, etc.) 
which are justified finally by the good agreement of the theory with experiment. One can guess, 
that such additional hypothesis are in fact the consequences of the dynamical laws. This idea 
was realized in the Boltzmann's approach, for which ergodicity of the Hamiltonian flows plays an 
important role. Verification of the latter for real systems is a difficult problem. In [1 ergodicity 
of a system of N hard spheres was proved. However, this proof can not be extended to oscillatory 
systems. In many papers [HI El H] the ergodicity problem is analyzed analytically and numerically. 

For quantum systems, justification of the micro-canonical distribution is even more problematic, 
since the matrix elements of the density operator calculated in the basis of the stationary states 
only change their phases and do not approach to the equilibrium distribution. Nevertheless, von 
Neumann [5] managed to prove some important statements, which could be treated as quantum 
basis of the classical statistical mechanics. In the von Neumann's approach, it is allowed to measure 
not all dynamical quantities, but only classical (macroscopic) variables which quantum operators 
mutually commute. 

One can guess, that in a non-isolated system being in contact with the thermal bath, all 
dynamical variables are characterized by statistical properties. In other words, the density operator 
of such a system at almost all time moments is the Gibbs distribution. In the present paper 
this assumption is verified for the linear chain of harmonic oscillators - the system of interacting 
particles, which allows exact solution of the Schrodinger equation. The dynamics of relaxation in 
this system is studied as well. 

Ukrainian Physical Journal, 25, N 7, 1135-1141 (1980). 
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1 Model 



We consider the periodical chain of harmonic oscillators. Each of them is coupled with two neigh- 
bours and with its equilibrium point |6]. The Hamiltonian has the form 

JV-l 

H = E 2 [pt + U(\qk-qk+l\) + n 2 Q q 2 k ] , q N = q . (A.l) 

k=0 

Here qk, pk are the coordinate and momentum of the fc-th oscillator, k = 0, . . . , N — 1; N is the 
numbers of oscillators, U(\qk — qk+i\)/2 is the interaction potential of the neighbouring oscillators, 
f2o is their own frequency, the Plank constant and the oscillator mass are put to the unit value, 
periodical boundary conditions are chosen. The fc-th oscillator will be called 'the fc-th atom'. 
Model (|A.1[) has been studied in [7] in order to calculate correlation functions and to estimate the 
Poincare periods in the thermodynamics limit. 

The dynamical problem can be solved exactly if U{\qk — <7fc+i|) = ^ 2 (<Zfc — <Zfc+i) 2 - I n this case 
the normal coordinate Q T , P T are introduced as 

q k = J2 eikT Qr/^N, p k = J2 e ~ ikTp -r/^> ( A - 2 ) 

T T 

t = 2itI/N, -N/2<l<N/2, 
with integer I. Two sets of secondary quantization operators are related with q,p, Q, P: 

t n q k + ipk , , qk - Wk 



y2T2g y2T2o 

_ UJ T Qt + iP\ f _ UrQj- iPr ,. , 

V 2w T V 2w T 

where = fig + ft 2 (2 sin §) 2 . 

The Hamiltonian can be rewritten in the form 

H = J^(4 a T + l/2)w r , (A.4) 

T 

which shows that the normal modes (phonons) do not interact. The eigenstate of the Hamiltonian 
reads as 

\m,...,n N ) = (ni\... njv!)- 1/2 (4r . . . (a^)"* |0), (A.5) 

where |0) is the ground state, is the number of phonons with the quasi-momentum Tj, Tj < r i+1 . 
If f2 <C f2o> each atom can be considered as a subsystem which interacts weakly with surrounding 
(the thermal bath). 



2 Relaxation of the one-particle density operator 

If fl is small, one can follow the relaxation dynamics of the single-atom density operator for the 
system prepared in some realistic non-equilibrium initial state. 
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Let us introduce the following notations; 

a n = JJflr T , n= {n T }, \n) = \n\ .. t ,n^), n! = JJn r !, 

r r 

77 I 

Ci-i^. etc. (A.6) 

If the whole system is described by the density operator p = p nn >\n){n'\, than the atom with 
the number k = in the chain has the density operator Tr^'p = p nn tTr^\n)(n'\, where 

CO CO oo 

Tr^A = E E , , ,, , ^r°l°> >< 

m=0m'=0m 1 ,..., row _ 1 =0 V m oW m l ! • • • m N~l<- 

(O^fo™ 1 ...A... (b\) mi {bl)< |0) (016™°. 

We have calculated the quantities Tr^ \n)(n'\ only at il = 0: 



.7 P £ T ^t- — P 

where |J) = ^^f-|0), 91 = £ T ™r = E r <■ The state I") has the 

energy E n = Eq + Slo^Tt, where 

Eo is the ground state energy. 

Using (IA.7I) it is possible to show that the atom's density operator which corresponds to the 
stationary state \n) of the system, has the form of the Gibbs distribution in the thermodynamic 
limit (N oo, Ot/iV = n = const) at least for n < 1, if the number aN of nonzero integers in the 
set n = {n T } is macroscopic (i.e. a ~ 1): 

Tr«|n)(n| = J— £ f_^) | J) ( J| [l + O^ 1 )] . (A.8) 
1 + n ' \ 1 + n J 

Perhaps, this result holds to some extent in other dynamical systems. 

Let us show, how the evolution of the system from some non-equilibrium initial state can be 
described. Let 

'0, t < 0, 

<2(t) = {CI, < t < T, (A.9) 
0, t>T, 

and f2 is so small, that one can neglect the terms of order e in 

(n'[fi(0 = n[Q(i) = 0] ) = 5 nn , + w, 

i.e. interaction almost does not change the stationary states \n). At t < 0, t > T one can use 
expression (|A.7[) . and at < t < T the system relaxes due to the weak interaction. 

If at t — the system is in the stationary state |n), then at time T its state vector is given by 
\n) exp[— i(ruj)T] (the terms of order e are omitted), where (am) = ^ T w r n T . 
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Let us choose the non-equilibrium initial state as 



(*>S 



1 



01! 



|V>(t)|*=o) = ^r=|0) = 



(A.10) 



the 0-th atom is exited to the energy OIJlo, while all the rest atoms remain unexcited. Then one 
gets from (|A~7l) . (fP0|) 



01! 



P! 2 (-l) 



P-J 



1 



E 5 



N p JUP~J)l ^ 11 



Eexp{i[fc(nr) — t(nuj)}} sr-^ 
(n - 01 ^ 



3xp{i[— fc(n'r) +t(n'ui)]} 
(n'-l)\ : 



(A.11) 



where p(k,t) is the density operator of the fc-th atom (A; = 0, 1, . . . , N — 1) at the time t, and 
(jit) — n T T. Straightforward calculations yield 



P (k,t) = j2\J)(J\ 



oi! 



J!(01 — J)! 



2.7 



(A.12) 



where 



KM) = 



exp i(fcr — ui T t) 



Keeping the linear terms in (il/ilo) 2 in the expansion of uj t , and replacing the sum in r in 
r(q, t) by the integral, we come to the integral representation of the Bessel function of the g-th 
order: 



KM) = \NJ k [n (n/n ) 2 t} 



(A.13) 



One can safely replace the sum by the integral at large N and if t <C A r /[57o(il/fio) 2 ] (the period 
of the integrand is much large than At = 2n/N). 

One can easily see, that the energy of the fc-th atom at the time t exceeds the ground state 
energy by the quantity 



(A.14) 



Perturbation reaches the fc-th atom after the time of order k/[fto (fl/fto)} 2 , and after that 
the energy of this atom oscillates slowly with decreasing amplitude. When the excitation energy 
distributes throughout the whole chain, one can not any more replace the sum (TA~T2l) by the 
integral. For the classical harmonic chain, formula (|A.14[) was obtained in [7]. 
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3 Time average and dispersion of one-particle density op- 
erator 



Let us introduce the following notations for the average quantities: 



(A) 



p(t) 



lim — 

T^oo T 



Tr Ap(t), 

T 

f(t)dt, 



where A, f are an operator a function of time, respectively. 

In order to calculate (A) p , it is sufficient to know p. Consider the space of phases with points 
ip = (tpi, . . . ,(ppf), < (fii < 2ir (we use the same notations for indexes as in formula (|A.5[) . and 
the measure dip = Y[ T dip T /(2ir). Let us define the operator-valued function p(k,ip) on this space 
which is obtained from (jA.111) by the replacement 



t(nuj) (nip) = ^""^ n T ip T , t(n'uj) (n'ip) 

T 

Due to the equalitj0 



n T ip T 



e^jp[it(u>(n' — n))} — (e~x.p[i(ip(n — n))} 



(A.Ua) 



one can conclude from (IA.11[) . that the time average of the operator p(k,t) equals to the phase 
average of the operator p(k, ip). This allows us to write the time average of the operator (|A.12[) in 
the form 



p = J2 \J)(J\ Cm {[r^)/N] 2J |l - [r(p)/N] 
J 



m-j 



(A.15) 



where 



r(ip) 



^2exp(iip T ) 



The problem is reduced to the random walk problem, which was studied in [S]. 
(51)-(53), (103), (104) of reference [5] we find 



Using formulas 



p N (r) = (S(r-r(cp))) lp =r / dppJ (p) J a (r p), 

Jo 

PN{r) -> — e r 1 , at N -> oo. 



(A.16) 



2 [S.R. 2012] The first equality in equation (A. 14a) is wrong, since it does not take into account that the energy 
spectrum E n = (u)n) of model jATTj is strongly degenerate. The corrected form of this equation is given by formula 
341) in Subsection 12.31 This correction induces also certain modifications into the calculation of p in this page and 
into the final result llA,17j l, These modifications are described in Proposition 12.21 and equation l|35|l in Subsection 
12.31 Fortunately, the above mentioned modifications are minimal, and equation JA.17I I holds for all but one or 
two atoms in the chain, depending on the parity of N . On the other hand, application of a weak anharmonic 
interaction breaking the energy spectrum degeneracy restores the result HA. 171 1: the time average of the reduced 
density operators of all atoms in the chain is given by the Gibbs distribution with the same temperature. 
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Using the above large- N asymptotics for pj\r(r) and for the asymptotic formula 

[l- (r/N) 2 ]* n ~ J ^ exp[-r 2 m/N 2 }, at N -> oo, 
we get 

lim p N = Y — L_ ( -^—] \J)(J\. (A.17) 

The canonical distribution is obtained for all values of n. 

It is clear from (|A.12[> that the time fluctuations of the density operator do not vanish in the 
thermodynamic limit. This behavior is caused by the degeneracy of the system in the differences 
of the energy levels of the stationary states: 

(un) - (urn 1 ) = (oj(ti + m)) - (w(n' + m)). (A. 18) 

If a weak enough anharmonic interaction is applied 

(a) this will not change substantially p, 

(b) at least at n < 1, the dispersion 



Djj>(p) = (pj - pj)*{pj- - pj>) 

vanishes in the thermodynamic limit. 

A weak anharmonic interaction can be described as a perturbation leading to a small shift of each 
energy level 

E n = (urn) + 6E n . 

Therefore, in equation (jA.llj) . we should replace now (nu) by E n , and (n'uj) by E n >, and the 
statement (a) becomes evident. 

Let us prove (b). One obtains from (| A. 12j) 

TO 2 p?(-l)P-r Q!2(_!)Q-J' 

JJ 'W N 2m 2^ n p J\(P - J)\ NQ J'\(Q - J')\ X 

E l\(n - l)l(n' - 1)1 k\(n - k)\(n' - k)V 

J2 T i T =m-p.j: T k r =<n-Q 

In the sum (|A.19|i . the following inequality holds \\n — n'\\ < 2min(P,Q), where ||n|| = 52 \ n r\- 
Since 

m\ 2 



E 



< 



N 2m ^ n(n-l)\M-l)\ kUn - k)Uri - k)\ 

c * c m N 2m n]n n ^ P \q\ ' (A.2U) 

0<jjn-n'jj<2min(P,Q) 



2G 



the series in (IA.19|) uniformly converges at n < 1. Therefore, one can treat P and Q as finite 
numbers. Let P < Q. The number of terms in the sum 



0<||n-ra'||<2P 

is smaller than 



y - 



(N + P -1)1 
P\(N-1)\ 



therefore 

ot! 2 



N 2m 



E 



E 



i w 2 
— < 



=01 



E T <=oi 

0<||n-n'||<2P 



N 2< ^ 



E 



=91 



(2V + P-1)! 
P!(JV- 1)! 



-1 2 



< 



p|2 • 



(A.21) 



which proves statement (b). It is quite possible that (IA.8I) and (b) hold also for n > 1. 

In conclusion let us note that the considered model can be useful for many other applica- 
tions. Since in the case of the harmonic interaction, its dynamics is integrable, the model can be 
used to verify different approximate methods of studying relaxation processes. In the case of the 
anharmonic interaction, one can study analytically or numerically the relaxation dynamics in a 
nonlinear system. The model could be useful also for analyzing the arising statistical properties in 
a non-macroscopic systems, like three-atom molecules. 

The author is grateful to S. V. Pclctminski for the discussion of the results of this work. 
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